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1. Introduction 

The CLEO Collaboration has recently reported the observation of some rare two- 
body decays of the type B — > irK, as well as interesting upper bounds for the decays 
B — > 7T7T and B — > KK [1]. In particular, they find the CP-averaged branching ratios 



Br(E -> 71- K + ) + Br(5 -> 7r + ^")l = (1.4 ± 0.3 ± 0.1) x 10~ 5 , 



Br(E + -> vr+fT ) + Br(£ _ -> tt"^ 



1.4 ±0.5 ±0.2) x 10~ 5 , 



Br(5 + -> 7r°ir + ) ± Br( J B~ -> vr ^)! = (1.5 ± 0.4 ± 0.3) x 10 



\-5 



(1.1) 
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This observation caused a lot of excitement, because these decays offer interesting 
insights into the relative strength of various contributions to the decay amplitudes, 
whose interference can lead to CP asymmetries in the decay rates. It indeed appears 
that there may be potentially large interference effects, depending on the magnitude 
of some strong interaction phases (see, e.g., [2]). Thus, although at present only 
measurements of CP-averaged branching ratios have been reported, the prospects 
are good for observing direct CP violation in some of the B — > ixK or B — > KK 
decay modes in the near future. 

It is fascinating that some information on CP-violating parameters can be ex- 
tracted even without observing a single CP asymmetry, from measurements of CP- 
averaged branching ratios alone. This information concerns the angle 7 of the 
so-called unitarity triangle, defined as 7 = &ig[(V* h V u d) / {V* b V c d)}. With the stan- 
dard phase conventions for the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 7 = 
arg(V^) to excellent accuracy. There have been proposals for deriving bounds on 7 
from measurements of the ratios 

t{B+) Br(B° -> 1T-K+) + Brf^B -> n+K~) 



R 



t(B°) Br(B+ -> n+K°) + Bt(B~ -> ir-R ) ' 



Br(E+ -> Tr+if °) + Bt(B- -> 71 K°) 
* 2[Br(5+ -> tt°K+) + Br(B- -> vr ^")] ' 1 ' J 

whose current experimental values are i? = 1.07 ± 0.45 (we use t(B + )/t(B°) = 
1.07 ± 0.03) and R* = 0.47 ± 0.24. The Fleischer-Mannel bound R > sin 2 7 [3] 
excludes values around |7| = 90° provided that R < 1. However, this bound is 
subject to theoretical uncertainties arising from electroweak penguin contributions 
and strong rescattering effects, which are difficult to quantify [4] -[9]. The bound 

1 - \[rZ < e 3 / 2 \S E w ~ cos 7I + 0(63/2) (1.3) 

derived by Rosner and the present author [10], where <5ew — 0.64 ±0.15 accounts for 
electroweak penguin contributions, is less affected by such uncertainties; however, it 
relies on an expansion in the small parameter 



£3/2 = y/2 Rsv(3) tan^ c 



Bl(B+ -> 7T + 7T ) + Bl(B- -> 7T-7T ) 1 ^ 



(1.4) 



Bi(B+ -> TT+iT ) + Br(B- -> tt-K ) 

whose value has been estimated to be £3/2 = 0.24±0.06. Here is the Cabibbo angle, 
and the factor i?su(3) — fx/ fn accounts for SU(3)-breaking corrections. Assuming 
the smallness of certain rescattering effects, higher-order terms in the expansion in 
£3/2 can be shown to strengthen the bound (1.3) provided that the value of i?* is not 
much larger than indicated by current data, i.e., if R* < (1 — e 3 / 2 /v / 2) 2 ~ 0.7 [10]. 

Our main goal in the present work is to address the question to what extent these 
bounds can be affected by hadronic uncertainties such as final-state rescattering ef- 
fects, and whether the theoretical assumptions underlying them are justified. To this 
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end, we perform a general analysis of the various B — > ttK decay modes, pointing 
out where theoretical information from isopsin and SU(3) flavour symmetries can be 
used to eliminate hadronic uncertainties. Our approach will be to vary parameters 
not constrained by theory (strong-interaction phases, in particular) within conserva- 
tive ranges so as to obtain a model-independent description of the decay amplitudes. 
An analysis pursuing a similar goal has recently been presented by Buras and Fleis- 
cher [11]. Where appropriate, we will point out the relations of our work with theirs 
and provide a translation of notations. We stress, however, that although we take a 
similar starting point, some of our conclusions will be rather different from the ones 
reached in their work. 

In Section 2, we present a general parametrization of the various isospin ampli- 
tudes relevant to B — > nK decays and discuss theoretical constraints resulting from 
flavour symmetries of the strong interactions and the structure of the low-energy 
effective weak Hamiltonian. We summarize model-independent results derived re- 
cently for the electroweak penguin contributions to the isovector part of the effective 
Hamiltonian [6, 9] and point out constraints on certain rescattering contributions 
resulting from B — > KK decays [7, 9, 12, 13]. The main results of this analysis 
are presented in Section 2.6, which contains numerical predictions for the various 
parameters entering our parametrization of the decay amplitudes. The remainder of 
the paper deals with phenomenological applications of these results. In Section 3, we 
discuss corrections to the Fleischer-Mannel bound resulting from final-state rescat- 
tering and electroweak penguin contributions. In Section 4, we show how to include 
rescattering effects to the bound (1.3) at higher orders in the expansion in £3/2- De- 
tailed predictions for the direct CP asymmetries in the various B — > nK decay modes 
are presented in Section 5, where we also present a prediction for the CP-averaged 
B° — > 7i°K° branching ratio, for which at present only an upper limit exists. In 
Section 6, we discuss how the weak phase 7, along with a strong-interaction phase 
difference 0, can be determined from measurements of the ratio -R* and of the di- 
rect CP asymmetries in the decays B ± — > -k^K^ and B ± — > ^K (here K° means 
K° or K°, as appropriate). This generalizes a method proposed in [14] to include 
rescattering corrections to the B^ — > n ± K° decay amplitudes. Section 7 contains a 
summary of our result and the conclusions. 

2. Isospin decomposition 
2.1 Preliminaries 

The effective weak Hamiltonian relevant to the decays B — > ttK is [15] 

n = -^\ £ C^XuQ? + XcQ*) - \ t J2CiQi +h.c, (2.1) 
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where X q = V* b V qs are products of CKM matrix elements, Cj are Wilson coefficients, 
and Qi are local four-quark operators. Relevant to our discussion are the isospin 
quantum numbers of these operators. The current-current operators Q\ 2 ~ bsuu 
have components with AI = and AI = 1; the current-current operators Q\ 2 ~ bscc 
and the QCD penguin operators Qz,...,e, ~ bsJ^QQ have AI = 0; the electroweak 
penguin operators <57,...,io ~ ^ s e g<?<L where e q are the electric charges of the quarks, 
have AI = and AI = 1. Since the initial i? meson has I = \ and the final states 
(7r.fr) can be decomposed into components with I — | and / = |, the physical 
— > 7ri^ decay amplitudes can be described in terms of three isospin amplitudes. 
They are called Bi/ 2 , Ay 2 , and A 3 / 2 referring, respectively, to AI = with I^k = \, 
AI = 1 with I^k = \, and AI = 1 with I wK = | [6, 16, 17]. The resulting expressions 
for the decay amplitudes are 



A(B + - 


-> n + K°) 


= -Bl/2 


+ A 1/2 


+ A 3/2 , 


V2A(B + - 


-> 7l°K + ) 


= ^1/2 


+ A 1/2 


- 2A 3/2 


-A(B° - 


> 7T-^ + ) 


= Bi/2 


-A 1/2 


- A 3 / 2 , 


V2A{B° - 


^ TT ^ ) 


= Bi/2 


-A 1/2 


+ 2A3/2 



From the isospin decomposition of the effective Hamiltonian it is obvious which 
operator matrix elements and weak phases enter the various isospin amplitudes. 
Experimental data as well as theoretical expectations indicate that the amplitude 
B1/2, which includes the contributions of the QCD penguin operators, is significantly 
larger than the amplitudes Ai/ 2 and A 3 / 2 [2, 6]. Yet, the fact that Ai/ 2 and A 3 / 2 are 
different from zero is responsible for the deviations of the ratios R and R* in (1.2) 
from 1. 

Because of the unitarity relation \ u + \ c + \ t = there are two independent CKM 
parameters entering the decay amplitudes, which we choose to be 1 — A c = e l7r |A c | and 
X u = e n |A u |. Each of the three isospin amplitudes receives contributions proportional 
to both weak phases. In total, there are thus five independent strong-interaction 
phase differences (an overall phase is irrelevant) and six independent real amplitudes, 
leaving as many as eleven hadronic parameters. Even perfect measurements of the 
eight branching ratios for the various B — > irK decay modes and their CP conjugates 
would not suffice to determine these parameters. Facing this problem, previous au- 
thors have often relied on some theoretical prejudice about the relative importance of 
various parameters. For instance, in the invariant SU(3)-amplitude approach based 
on flavour-flow topologies [18, 19], the isospin amplitudes are expressed as linear com- 
binations of a QCD penguin amplitude P, a tree amplitude T, a colour-suppressed 
tree amplitude C, an annihilation amplitide A, an electroweak penguin amplitude 

1 Taking A c to be real is an excellent approximation. 
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Pew; an d a colour-suppressed electroweak penguin amplitude P^ w , which are ex- 
pected to obey the following hierarchy: |P| 3> \T\ ~ |-Pew| ^> |C| ~ |PewI > 1^1- 
These naive expectations could be upset, however, if strong final-state rescattering 
effects would turn out to be important [5]-[8], a possibility which at present is still 
under debate. Whereas the colour-transparency argument [20] suggests that final- 
state interactions are small in B decays into a pair of light mesons, the opposite 
behaviour is exhibited in a model based on Regge phenomenology [21]. For compar- 
ison, we note that in the decays B — > D^h, with h = it or p, the final-state phase 
differences between the I — \ and I — | isospin amplitudes are found to be smaller 
than 30°-50° [22]. 

Here we follow a different strategy, making maximal use of theoretical constraints 
derived using flavour symmetries and the knowledge of the effective weak Hamilto- 
nian in the Standard Model. These constraints help simplifying the isospin amplitude 
A3/2, for which the two contributions with different weak phases turn out to have 
the same strong-interaction phase (to an excellent approximation) and magnitudes 
that can be determined without encountering large hadronic uncertainties [10]. The- 
oretical uncertainties enter only at the level of SU(3)-breaking corrections, which can 
be accounted for using the generalized factorization approximation [22]. Effectively, 
these simplifications remove three parameters (one phase and two magnitudes) from 
the list of unknown hadronic quantities. There is at present no other clean theo- 
retical information about the remaining parameters, although some constraints can 
be derived using measurements of the branching ratios for the decays B^ — > K^K 
and invoking SU(3) symmetry [7, 9, 12, 13]. Nevertheless, interesting insights can 
be gained by fully exploiting the available information on A 3 / 2 . 

Before discussing this in more detail, it is instructive to introduce certain linear 
combinations of the isospin amplitudes, which we define as 

B 1/2 + A 1/2 + A 3/2 = P + A- -P° w , 

—3A 3 / 2 = T + C + Pew + Pew > 
-2(A 1/2 + A 3/2 ) = T-A + Pg w . (2.3) 

In the latter two relations, the amplitudes T, C and A carry the weak phase e n , 
whereas the electroweak penguin amplitudes Pew an d Pew carr y the weak phase 2 
e l7T . Decomposing the QCD penguin amplitude as P = J2 q \Pq, an d similarly writing 
A = \ U A U and P^ w = \tP%^,ti we rewrite the first relation in the form 

P1/2 + A 1/2 + A 3/2 = -A C (P - P c - §P|w it ) + \u(A u -P t + P u ) 

2 Because of their smallness, it is a safe approximation to set X t = — X c for the electroweak 
penguin contributions, and to neglect electroweak penguin contractions in the matrix elements of 
the four-quark operators Qf and Q\. 
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= \P\e i<t,p (e in + e a e il e ir >y (2.4) 

By definition, the term \P\ e"^ p e l7r contains all contributions to the B + — > ir + K° de- 
cay amplitude not proportional to the weak phase e n . We will return to a discussion 
of the remaining terms below. It is convenient to adopt a parametrization of the 
other two amplitude combinations in (2.3) in units of \P\, so that this parameter 
cancels in predictions for ratios of branching ratios. We define 



-3^3/2 
|P| 

-2{A 1/2 + A 3/2 ) 



e T e i4>T (e 11 -q c e itJJC ) , (2.5) 



where the terms with q and qc arise from electroweak penguin contributions. In the 
above relations, the parameters rj, ^3/2, 4>t, w, and ojc are strong-interaction phases. 
For the benefit of the reader, it may be convenient to relate our definitions in (2.4) 
and (2.5) with those adopted by Buras and Fleischer [11]. The identificantions are: 
\P\ e t<i>p <-> A c |P te | e l5tc , e a e tv —pe te , 3 / 2 <-> St+c, and 4>t <-> St- The notations 
for the electroweak penguin contributions conincide. Moreover, if we define 

e-3/2 = t = £V2 , (2.6) 

J 1 — 2e a cos 77 cos 7 + 

then e 3 / 2 <-> r c and erf £3/2 ^ r/r c . With this definition, the parameter e 3 / 2 is 
precisely the quantity that can be determined experimentally using the relation (1.4). 

2.2 Isovector part of the effective weak Hamiltonian 

The two amplitude combinations in (2.5) involve isospin amplitudes defined in terms 
of the strong-interaction matrix elements of the AI = 1 part of the effective weak 
Hamiltonian. 3 This part contains current-current as well as electroweak penguin op- 
erators. A trivial but relevant observation is that the electroweak penguin operators 
Qg and Q10, whose Wilson coefficients are enhanced by the large mass of the top 
quark, are Fierz-equivalent to the current-current operators Qi and Q2 [6, 9, 23]. As 
a result, the AI = 1 part of the effective weak Hamiltonian for B — > irK decays can 
be written as 

n A i=i = ^ { (a u Ci - ^XtC^j Qi + (\ U C2 - ^A t C 10 ) Q 2 + ■ ■ j + h.c. , (2.7) 

where Qi = \{Q^ — Qf) are isovector combinations of four-quark operators. The dots 
represent the contributions from the electroweak penguin operators Qj and Q$, which 



3 This statement implies that QED corrections to the matrix elements are neglected, which is an 
excellent approximation. 
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have a different Dirac structure. In the Standard Model, the Wilson coefficients of 
these operators are so small that their contributions can be safely neglected. It is 
important in this context that for heavy mesons the matrix elements of four-quark 
operators with Dirac structure (V — A) <g> (V + A) are not enhanced with respect 
to those of operators with the usual (V — A) <g> (V — A) structure. To an excellent 
approximation, the net effect of electroweak penguin contributions to the AJ = 1 
isospin amplitudes in B — > nK decays thus consists of the replacements of the Wilson 
coefficients C\ and C 2 of the current-current operators with the combinations shown 
in (2.7). Introducing the linear combinations C± = (C 2 ±Ci) and Q± = ^(Q 2 ±Qi), 
which have the advantage of being renormalized multiplicatively, we obtain 

Hai=i * ^\V: b V us \ {c + (e* - 5 + ) Q + + C-(e* - 5-) <?_} + h.c. , (2.8) 

where 



3 cot 9 C C10 ± C 9 



2\V ub /V cb \ C 2 ±d 



(2.9) 



We have used \ u /\ t ^ — A u /A c ~ — tan 9c \V ub /V cb \ e' n , with the ratio [Kfc/Kbl = 
0.089 ± 0.015 determined from semileptonic B decays [24]. 

From the fact that the products C±Q± are renormalization-group invariant, it 
follows that the quantities 5± themselves must be scheme- and scale-independent (in 
a certain approximation). Indeed, the ratios of Wilson coefficients entering in (2.9) 
are, to a good approximation, independent of the choice of the renormalization scale. 
Taking the values d = -0.308, C 2 = 1.144, C 9 = -1.280a and C w = 0.328a, 
which correspond to the leading-order coefficients at the scale fi = m b [15], we 
find (C10 + C 9 )/(C 2 + Ci) w -1.14a and (C 10 - C 9 )/(C 2 - Ci) w 1.11a, imply- 
ing that 5- ~ —5 + to a good approximation. The statement of the approximate 
renormalization-group invariance of the ratios 6± can be made more precise by not- 
ing that the large values of the Wilson coefficients C 9 and C 10 at the scale // = m b 
predominantly result from large matching contributions to the coefficient Cg(mw) 
arising from box and Z-penguin diagrams, whereas the O(a) contributions to the 
anomalous dimension matrix governing the mixing of the local operators Qi lead 
to very small effects. If these are neglected, then to next-to-leading order in the 
QCD evolution the coefficients (C10 ± Cg) are renormalized multiplicatively and in 
precisely the same way as the coefficients {C 2 ± C\). We have derived this result us- 
ing the explicit expressions for the anomalous dimension matrices compiled in [15] . 4 
Hence, in this approximation the ratios of coefficients entering the quantities 5± are 



4 The equivalence of the anomalous dimensions at next-to-leading order is nontrivial because 
the operators Qg and Q10 are related to Q\ and Q2 by Fierz identities, which are valid only in 
four dimensions. The corresponding two- loop anomalous dimensions are identical in the naive 
dimensional rcgularization scheme with anticommuting 75 . 
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renormalization-scale independent and can be evaluated at the scale m-w, so that 
C W ±C 9 , „ , \ a x t f 3kix t \ 

n ~T~n ~ ±C 9 (m w ) = Tj^ . 2Q 1 + r + . . . ~ Tl-18a , (2.10) 

C 2 ± Ci 127T sin ^ y x t — 1 y 

where #vk is the Weinberg angle, and x t = {m t /m w ) 2 . This result agrees with 
an equivalent expression derived by Fleischer [23]. The dots in (2.10) represent 
renormalization-scheme dependent terms, which are not enhanced by the factor 
l/sin 2 #w These terms are numerically very small and of the same order as the 
coefficients C 7 and C 8 , whose values have been neglected in our derivation. The 
leading terms given above are precisely the ones that must be kept to get a consis- 
tent, renormalization-group invariant result. We thus obtain 

a cot On Xt ( 31nrr+\ . 

5 + = -6- = — t- t ,t , — t— 1 + = 0.68 ± 0.11 , (2.11) 

+ 8tt \V ub /V cb \ sm 2 9 w \ x t -lj ' 1 ; 

where we have taken a = 1/129 for the electromagnetic coupling renormalized at the 
scale m b) and m t = m t (m t ) = 170 GeV for the running top-quark mass in the MS 
renormalization scheme. Assuming that there are no large 0(a s ) corrections with 
this choice, the main uncertainty in the estimate of 5 + in the Standard Model results 
from the present error on \V u b\, which is likely to be reduced in the near future. 

We stress that the sensitivity of the B — > nK decay amplitudes to the value of 
5 + provides a window to New Physics, which could alter the value of this parameter 
significantly. A generic example are extensions of the Standard Model with new 
charged Higgs bosons such as supersymmetry, for which there are additional matching 
contributions to C${mw)- We will come back to this point in Section 4. 

2.3 Structure of the isospin amplitude A 3 / 2 

[/-spin invariance of the strong interactions, which is a subgroup of flavour SU(3) 
symmetry corresponding to transformations exchanging d and s quarks, implies that 
the isospin amplitude A 3 / 2 receives a contribution only from the operator Q + in 
(2.8), but not from Q_ [10]. In order to investigate the corrections to this limit, we 
parametrize the matrix elements of the local operators C±Q± between a B meson 
and the (irK) isospin state with I — | by hadronic parameters Ky 2 e^ 3 / 2 , so that 

-3A 3/2 = Kt /2 e</* ( e * - S+ ) + K~ /2 e»v> (e* + S + ) 

ee (K^e l ^ + K- /2 e^)(e^-qen. (2.12) 

In the SU(3) limit K^ 2 = 0, and hence SU(3)-breaking corrections can be parametrized 
by the quantity 



Ke i^ 3 /2 = _ 3/2 



^3/2 



(2.13) 
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in terms of which 

qd lw = (l-/te lA ^/ 2 )5 + . (2.14) 

This relation generalizes an approximate result derived in [10]. 

The magnitude of the SU(3)-breaking effects can be estimated by using the 
generalized factorization hypothesis to calculate the matrix elements of the current- 
current operators [22]. This gives 



k ~ 2 



ai + a 2 A K + A w ^ 



(6 ±6)%, A^ 3 / 2 ~0, (2.15) 

_ai - a 2 A K - A n J 

where Ax = fK(m B —m 2 r )F^^' n (mj < ) and A n = f w (m B — m 2 K )F^^ K (m 2 r ) are combi- 
nations of hadronic matrix elements, and a\ and a 2 are phenomenological parameters 
defined such that they contain the leading corrections to naive factorization. For a 
numerical estimate we take a 2 /a 1 = 0.21 ±0.05 as determined from a global analysis 
of nonleptonic two-body decays of B mesons [22], and A T /A K = 0.9 ± 0.1, which is 
consistent with form factor models (see, e.g., [25]-[27]) as well as the most recent 
predictions obtained using light-cone QCD sum rules [28]. Despite the fact that 
nonfactorizable corrections are not fully controlled theoretically, the estimate (2.15) 
suggests that the SU(3)-breaking corrections in (2.14) are small. More importantly, 
such effects cannot induce a sizable strong-interaction phase uj. Since Q + and Q_ 
are local operators whose matrix elements are taken between the same isospin eigen- 
states, it is very unlikely that the strong-interaction phases (f>^, 2 and <p^ 2 could differ 
by a large amount. If we assume that these phases differ by at most 20°, and that 
the magnitude of k is as large as 12% (corresponding to twice the central value 
obtained using factorization), we find that \u\ < 2.7°. Even for a phase difference 
A923/2 — \<f>2/2 ~ 03/2I = 90°, which seems totally unrealistic, the phase \u\ would not 
exceed 7°. It is therefore a safe approximation to work with the real value [10] 

S EW = (1 - k)5 + = 0.64 ±0.15, (2.16) 

where to be conservative we have added linearly the uncertainties in the values of 
k and S + . We believe the error quoted above is large enough to cover possible 
small contributions from a nonzero phase difference A</? 3 / 2 or deviations from the 
factorization approximation. For completeness, we note that our general results 
for the structure of the electroweak penguin contributions to the isospin amplitude 
A3/2, including the pattern of SU(3)-breaking effects, are in full accord with model 
estimates by Deshpande and He [29]. Generalizations of our results to the case of 
B — > 7T7T, KK decays and the corresponding B s decays are possible using SU(3) 
symmetry, as discussed in [30, 31]. 

In the last step, we define K+ /2 + j£- ^3/2 = |p| £3/2 e ^ 3/ 2 j so that [10] 

4^ = £ 3/2 e»*/*(e* - few) ■ (2.17) 
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The complex quantity qe lul in our general parametrization in (2.5) is now replaced 
with the real parameter <5ew> whose numerical value is known with reasonable accu- 
racy. The fact that the strong-interaction phase oo can be neglected was overlooked 
by Buras and Fleischer, who considered values as large as \uo\ = 45° and therefore 
asigned a larger hadronic uncertainty to the isospin amplitude A 3 / 2 [11]. 

In the SU(3) limit, the product \P\ e 3 / 2 is determined by the decay amplitude 
for the process B ± — > ^-n through the relation 

\P\ £3/2 = V2 ^^tan^ c \A(B ± - 7r ± 7 r°) ] , (2.18) 



where 5 

Rew 



A^Siwt-" (2.19) 



is a tiny correction arising from the very small electroweak penguin contributions to 
the decays B ± -> t^tt . Here A = sin0 c « 0.22 and R t = [(1 - pf + rf} 1 ' 2 ~ 1 
are Wolfenstein parameters, and a is another angle of the unitarity triangle, whose 
preferred value is close to 90° [32]. It follows that the deviation of -Rew from 1 is of 
order 1-2%, and it is thus a safe approximation to set -Rew — 1- More important 
are SU(3)-breaking corrections, which can be included in (2.18) in the factorization 
approximation, leading to 

p _ fl i Ik «2 F^ K {ml) f K _ 
Ol + «2 Jtt cti + a 2 -r (m^J /„■ 
where we have neglected a tiny difference in the phase space for the two decays. 
Relation (2.18) can be used to determine the parameter e 3 / 2 introduced in (2.6), 
which coincides with e 3 / 2 up to terms of 0(e a ). To this end, we note that the CP- 
averaged branching ratio for the decays B^ 1 — > n^K is given by 



1 



By(B ± -> ^K ) = ^ [Bi(B + -> 7i + K°) + Br(£ _ -> 7r"/iT ) 

= |P| 2 (l - 2e a cosr/cos7 + ^) . (2.21) 

Combining this result with (2.18) we obtain relation (1.4), which expresses e 3 / 2 in 
terms of CP-averaged branching ratios. Using preliminary data reported by the 
CLEO Collaboration [1] combined with some theoretical guidance based on factor- 
ization, one finds e 3/2 = 0.24 ± 0.06 [10]. 

To summarize, besides the parameter Sew controlling electroweak penguin con- 
tributions also the normalization of the amplitude A 3 / 2 is known from theory, albeit 
with some uncertainty related to nonfactorizable SU(3)-breaking effects. The only re- 
maining unknown hadronic parameter in (2.17) is the strong- interaction phase (f) 3 / 2 . 
The various constraints on the structure of the isospin amplitude A 3 / 2 discussed 
here constitute the main theoretical simplification of B — > nK decays, i.e., the only 
simplification rooted on first principles of QCD. 

5 We disagree with the result for this correction presented in [11]. 
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2.4 Structure of the amplitude combination -Bi/ 2 + ^1/2 + -A3/2 

The above result for the isospin amplitude A 3 / 2 helps understanding better the struc- 
ture of the sum of amplitudes introduced in (2.4). To this end, we introduce the 
following exact parametrization: 



B 1/2 + A 1/2 + A, 



1/2 



I3/2 



e m e trpp 



(2.22) 



where we have made explicit the contribution proportional to the weak phase e* 7 
contained in A 3 / 2 . From a comparison with the parametrization in (2.4) it follows 
that 



where 



>3/2 



e a e" 
and A 



e^(V A -l) , 



(2.23) 



>l/2 



>3/2- 



Of course, this is just a simple 



reparametrization. However, the intuitive expectation that e a is small, because this 
terms receives contributions only from the penguin (P u — P t ) and from annihilation 
topologies, now becomes equivalent to saying that £ e tA is close to 1, so as to allow for 
a cancelation between the contributions corresponding to final-state isospin I — \ 
and / = I in (2.23). But this can only happen if there are no sizable final-state 
interactions. The limit of elastic final-state interactions can be recovered from (2.23) 
by setting £ = 1, in which case we reproduce results derived previously in [5, 6]. Be- 
cause of the large energy release in B — > nK decays, however, one expects inelastic 
rescattering contributions to be important as well [7, 21]. They would lead to a value 

From (2.23) it follows that 



^3/2 



1 - 2£ cos A + £ 2 



2v^ 



■ £3/2* 



1 



IV 



• 2 A 
sm- 



(2.24) 



where without loss of generality we define e a to be positive. Clearly, e a £3/2 
provided the phase difference A is small and the parameter £ close to 1. There are 
good physics reasons to believe that both of these requirements may be satisfied. 
In the rest frame of the B meson, the two light particles produced in B — > nK 
decays have large energies and opposite momenta. Hence, by the colour-transparency 
argument [20] their final-state interactions are expected to be suppressed unless there 
are close-by resonances, such as charm-anticharm intermediate states (DD S , J/ip K, 
etc.). However, these contributions could only result from the charm penguin [33, 34] 
and are thus included in the term \P\ e^ p in (2.4). As a consequence, the phase 
difference <f> = <p 3 / 2 — (ftp could quite conceivably be sizable. On the other hand, the 
strong phases 03/2 and (fii/ 2 in (2.22) refer to the matrix elements of local four-quark 
operators of the type bsuu and differ only in the isospin of the final state. We believe 
it is realistic to assume that |A| = \4>\/ 2 — 4>z/ 2 \ < 45°. Likewise, if the parameter £ 
were very different from 1 this would correspond to a gross failure of the generalized 
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factorization hypothesis (even in decays into isospin eigenstates), which works so 
well in the global analysis of hadronic two-body decays of B mesons [22]. In view 
of this empirical fact, we think it is reasonable to assume that 0.5 < £ < 1.5. With 
this set of parameters, we find that e a < 0.35£ 3 / 2 < 0.1. Thus, we expect that the 
rescattering effects parametrized by e a are rather small. 

A constraint on the parameter e a can be derived assuming [/-spin invariance of 
the strong interactions, which relates the decay amplitudes for the processes — > 
7r ± i^° and B ± -> ^K up to the substitution [7, 9, 13] 

V* b V ud ~ y , Ae-O^-AAe, (2.25) 

where A ~ 0.22 is the Wolfenstein parameter. Neglecting SU(3)-breaking corrections, 
the CP-averaged branching ratio for the decays B ± — > K ± K° is then given by 

Br(B ± -> K^K ) = ^ [Bt(B+ -> K + K°) + Br(£ _ -> K~K°) 

= \P\ 2 [A 2 + 2e a cos 77 cos 7 + (e a /A) 2 ] , (2.26) 

which should be compared with the corresponding result for the decays B ± — > ir^K 
given in (2.21). The enhancement (suppression) of the subleading (leading) terms 
by powers of A implies potentially large rescattering effects and a large direct CP 
asymmetry in B^ — > K ± K° decays. In particular, comparing the expressions for the 
direct CP asymmetries, 

Br(B + -> 7T+ K°) - Br(B- -> Tr^i? ) 2e a sin 77 sin 7 



A CP (7T + X L 



Br(5+ -> vr+fT ) + Br (5- -> vr-fsT ) 1 - 2e a cos 77 cos 7 + ^ 



4 ^+^on 2e a sin 77 sin 7 

AcpI^ a J = - x2 1 o 71 — 7v~> \ 2 - 27 ) 

X 2 + 2e a cos 77 cos 7 + (e a /A) 2 

one obtains the simple relation [9] 

AcAK + K«) = Br(B± -» tt^AT ) 

A C p(7r+ir°) Br(£± -> ^±K°) ' 1 ' ' 

In the future, precise measurements of the branching ratio and CP asymmetry in 
B^ — > K ± K° decays may thus provide valuable information about the role of rescat- 
tering contributions in B ± — > ti ± K decays. In particular, upper and lower bounds 
on the parameter e a can be derived from a measurement of the ratio 

Br(B* -> K ± K°) = A 2 + 2 gQ cosr / cos7 + (e Q /A) 2 
* Br(£± -> tt^ ) l-2£ a cos77cos 7 + £ 2 ! ' 1 ' 

Using the fact that Rk is minimized (maximized) by setting cos 77 cos 7 = —1 (+1), 
we find that 

A(v^-A)^^A(^+A) 



1 + X^Rk - I- \^R K 
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This generalizes a relation derived in [7]. Using data reported by the CLEO Col- 
laboration [1], one can derive the upper bound R K < 0.7 (at 90% CL) implying 
e a < 0.28, which is not yet a very powerful constraint. However, a measurement of 
the branching ratio for B ± — > K ± K° could improve the situation significantly. For 
the purpose of illustration, we note that from the preliminary results quoted for the 
observed event rates one may deduce the "best fit" value Rk ~ 0.15 (with very large 
errors!). Taking this value literally would give the allowed range 0.03 < e a < 0.14. 

Based on a detailed analysis of individual rescattering contributions, Gronau and 
Rosner have argued that one expects a similar pattern of final-state interactions in 
the decays B ± -> K±K Q and B° -> K^K^ [12]. One could then use the tighter 
experimental bound Br(B° — > K ± K T ) < 2 x 10~ 6 to obtain e a < 0.16. However, this 
is not a model- independent result, because the decay amplitudes for B° — > K ± K T 
are not related to those for B ± — > n^K by any symmetry of the strong interactions. 
Nevertheless, this observation may be considered a qualitative argument in favour of 
a small value of e a . 

2.5 Structure of the amplitude combination A 1 / 2 + A 3 / 2 

None of the simplifications we found for the isospin amplitude A 3 / 2 persist for the 
amplitude Ai/ 2 . Therefore, the sum A 1 / 2 + A 3 / 2 suffers from larger hadronic uncer- 
tainties than the amplitude A 3 / 2 alone. Nevertheless, it is instructive to study the 
structure of this combination in more detail. In analogy with (2.12), we parametrize 
the matrix elements of the local operators C±Q± between a B meson and the (ttK) 
isospin state with I — \ by hadronic parameters K^ 2 e^ 1 / 2 , so that 

-3A l/2 = K+ 2 e^fye* - 6+) + K~ /2 e^e* + 5+) . (2.31) 
Next, we define parameters e' and r by 

Y(l±^^K 2 + ^/ 2 ) • (2-32) 

This general definition is motivated by the factorization approximation, which pre- 
dicts that r ~ a 2 /ai = 0.21 ±0.05 is the phenomenological colour-suppression factor 
[22], and 

~ aiAK = o.84 ± 0.04 . (2.33) 
£3/2 aiA K + a 2 A n 

With the help of these definitions, we obtain 
-2(A 1/2 + A 3/2 ) e' r 



2 
+ 



(1 + r) e^/2( e ^ - 5 + ) + (1 - r) e 1 ^ (e^ + 5+) 
^ (>3/2 - e ^v») (e^ - S + ) , (2.34) 
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where we have neglected some small, SU(3)-breaking corrections to the second term. 
Nevertheless, the above relation can be considered a general parametrization of the 
sum A1/2 + A 3 / 2 , since it still contains two undetermined phases and magnitudes 
e' and r. 

With the explicit result (2.34) at hand, it is a simple exercise to derive expressions 
for the quantities entering the parametrization in (2.5). We find 



£t 6 t = -e l<p ^ 



y A ^/2 + l) +r (e lA ^/2_l) 



_l_ 2g3 / 2 f e ^3/2 _ e % 2 



tu>c 



> + 



r(e iA01 /2 + 1) + ( e iA ^/2 - 1) + 


4^3/2 

3e' 


e i(4> 3/2 -4>+ /2 ) _ j 


( e iA^ 1/2 + 1) + r ( e *A0 1/2 _ 1) + 


4^3/2 

3e' 


e i{4>3/2-4>t/i) _ 1 



(2.35) 



where A0!/ 2 = 0^/2 — ^1/2 • This result, although rather complicated, exhibits in 
a transparent way the structure of possible rescattering effects. In particular, it 
is evident that the assumption of "colour suppression" of the electroweak penguin 
contribution, i.e., the statement that qc = 0{r) [3, 9, 11, 35], relies on the smallness of 
the strong-interaction phase differences between the various terms. More specifically, 
this assumption would only be justified if 



|A0i /2 | < 2r = 25 c 



?3/2 



1/2 



< 



3r e' 



£3/2 



= 15°. 



(2.36) 



We believe that, whereas the first relation may be a reasonable working hypothesis, 
the second one constitues a strong constraint on the strong-interaction phases, which 
cannot be justified in a model-independent way. As a simple but not unrealistic model 
we may thus consider the approximate relations obtained by setting A0i/ 2 = , 
which have been derived previously in [6]: 

e T e i4,T ~ e'e^v* + ^ (e* 4 *'* - e^/^j , 



IU>C 



r + 


2^3/2 


e i{4>3/2-4>t /2 ) 


- 1 


3e' 






1 + 


2^3/2 


e i(4> 3 /2-4>i /2 ) 


- 1 


3e' 







(2.37) 



The fact that in the case of a sizable phase difference between the I = \ and I = | 
isospin amplitudes the electroweak penguin contribution may no longer be as small 
as 0(r) has been stressed in [6] but was overlooked in [9, 11]. Likewise, there is some 
uncertainty in the value of the parameter et, which in the topological amplitude 
approach corresponds to the ratio |T — A|/]P| [19]. Unlike the parameter £3/2, 
the quantities e' and r cannot be determined experimentally using SU(3) symmetry 
relations. But even if we assume that the factorization result (2.33) is valid and 
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take £3/2 = 0.24 and e' = 0.20 as fixed, we still obtain 0.12 < e T < 0.20 depending 
on the value of the phase difference (03/2 — 4>t/2)- Note that from the approximate 
expression (2.37) it follows that Et < e' provided that e'/e 3 / 2 > 2/3, as indicated by 
the factorization result. This observation may explain why previous authors find the 
value e' = 0.15±0.05 [36], which tends to be somewhat smaller than the factorization 
prediction e' « 0.20. 

2.6 Numerical results 

Before turning to phenomenological applications of our results in the next sec- 
tion, it is instructive to consider some numerical results obtained using the above 
parametrizations. Since our main concern in this paper is to study rescattering ef- 
fects, we will keep e 3 / 2 = 0.24 fixed and assume that e' /s 3 / 2 = 0.84 ± 0.04 and 
r = 0.21 ± 0.05 as predicted by factorization. Also, we shall use the factorization 
result for the parameter k in (2.15). 
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Figure 1: Real and imaginary parts of the quantities x = e a e l<yri ^ (left) and y = 
£t e ^T-^/2) (right) for different choices of hadronic parameters. 

For the strong-interaction phases we consider two sets of parameter choices: one 
which we believe is realistic and one which we think is very conservative. For the 
realistic set, we require that 0.5 < f < 1.5, 103/2 - 1 } 2 ) | < 45°, and |0 + - <f>j\ < 20° 
(with I — | or |). For the conservative set, we increase these ranges to < £ < 2, 
103/2 - 0$! < 90°, and |0+ - 07| < 45°. In our opinion, values outside these ranges 
are quite inconceivable. Note that, for the moment, no assumption is made about 
the relative strong-interaction phases of tree and penguin ampltiudes. We choose the 
various parameters randomly inside the allowed intervals and present the results for 
the quantities e a e lT) in units of e 1 ^, Et^ t in units of e^ 3 / 2 , and qtc) e luJ( -°) in units 
of S + in the form of scatter plots in Figures 1 and 2. The black and the gray points 
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correspond to the realistic and to the conservative parameter sets, respectively. The 
same colour coding will be used throughout this work. 

The left-hand plot in Figure 1 shows 
that the parameter e a generally takes rather 
small values. For the realistic parameter 
set we find e a < 0.08, whereas values up 
tp 0.15 are possible for the conservative 
set. There is no strong correlation be- 
tween the strong-interaction phases i] and 
0. An important implication of these ob- 
servations is that, in general, there will 
be a very small difference between the 
quantities e 3 / 2 and £ 3 / 2 in (2.6). We shall 
therefore consider the same range of val- 
ues for the two parameters. From the 
right-hand plot we observe that for re- 
alistic parameter choices 0.15 < et < 
0.22; however, values between 0.08 and 
0.24 are possible for the conservative pa- 
rameter set. Note that there is a rather 
strong correlation between the strong-intera^/c^f) e 
phases 0t and 03/2, which differ by less 
than 20° for the realistic parameter set. We will see in Section 5 that this implies 
a strong correlation between the direct CP asymmetries in the decays B ± — > -n^K^ 
and B° — > it T K ± . Figure 2 shows that, even for the realistic parameter set, the 
ratio qc/8+ can be substantially larger than the naive expectation of about 0.2. In- 
deed, values as large as 0.7 are possible, and for the conservative set the wide range 
< qc/&+ < 1-4 is allowed. Likewise, the strong-interaction phase ojc can naturally 
be large and take values of up to 75° even for the realistic parameter set. (Note that, 
without loss of generality, only points with positive values of ui(c) are displayed in the 
plot. The distribution is invariant under a change of the sign of the strong-interaction 
phase.) This is in stark contrast to the case of the quantity qe luJ entering the isospin 
amplitude A 3 / 2 , where both the magnitude q and the phase u are determined within 
very small uncertainties, as is evident from the figure. 



Figure 2: Real and imaginary parts of the 
quantity z = (qc/5+) e tuJC . The accumula- 
tion of points in the lower-right corner shows 
the corresponding results for the quantity 
tuJ . Only points with Imz > are 

shown. 



3. Hadronic uncertainties in the Fleischer— Mannel bound 

As a first phenomenological application of the results of the previous section, we 
investigate the effects of rescattering and electroweak penguin contributions on the 
Fleischer-Mannel bound on 7 derived from the ratio R defined in (1.2). In general, 
R 7^ 1 because the parameter et in (2.5) does not vanish. To leading order in the 
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small quantities £j, we find 



R ~ 1 — 2e T cos cos 7 — q c cos((f) + wc) + 0(e 2 ) , 



(3.1) 



where = 0t — 0p. Because of the uncertainty in the values of the hadronic 
parameters et, qc and u>c, it is difficult to convert this result into a constraint 
on 7. Fleischer and Mannel have therefore suggested to derive a lower bound on the 
ratio R by eliminating the parameter Et from the exact expression for R. In the limit 
where E a and qc are set to zero, this yields R > sin 2 7 [3] . However, this simple result 
must be corrected in the presence of rescattering effects and electroweak penguin 
contributions. The generalization is [9] 



R > 



1 - 2q c E a cos(^c* + V) + Qc £ l 



(1 — 2q c cosuuc cos 7 + q c ){l — 2E a cos rj cos 7 + e%) 



sin 2 7 . 



(3.2) 
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The most dangerous rescattering effects arise from the terms involving the elec- 
troweak penguin parameter qc- As seen from Figure 2, even restricting ourselves to 
the realistic parameter set we can have 2qc cosuc ~ 5+ ~ 0.7 and q c ~ 0.5 5\ ~ 0.2, 
implying that the quadratic term in the denominator by itself can give a 20% cor- 
rection. The rescattering effects parametrized by £ a are presumably less important. 

The results of the numerical analy- 
sis are shown in Figure 3. In addition to 
the parameter choices described in Sec- 
tion 2.6, we vary e 3 / 2 and 5 + in the ranges 
0.24 ± 0.06 and 0.68 ± 0.11, respectively. 
Now also the relative strong-interaction 
phase between the penguin and / = 
I tree amplitudes enters. We allow val- 
ues |0| < 90° for the realistic parameter 
set, and impose no constraint on at all 
for the conservative parameter set. The 
figure shows that the corrections to the 
Fleischer-Mannel bound are not as large 
as suggested by the result (3.2), the rea- 
son being that this result is derived allow- 
ing arbitrary values of et, whereas in our 
analysis the allowed values for this pa- 
rameter are constrained. However, there 
are sizable violations of the naive bound 
R < sin 2 7 for I7I in the range between 




75. 100. 125. 150. 175. 
gamma 



Figure 3: Results for the ratio R versus I7I 
(in degrees) for different choices of hadronic 
parameters. The curve shows the minimal 
value -Rmm = sin 2 7 corresponding to the 
Fleischer-Mannel bound. The band shows 
the current experimental value of R with its 
la variation. 



65° and 125°, which includes most of the 
region 47° < 7 < 105° preferred by the global analysis of the unitarity triangle [32]. 
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Whereas these violations are numerically small for the realistic parameter set, they 
can become large for the conservative set, because then a large value of the phase dif- 
ference |03/2 — 0jy 2 1 1S allowed [6]. We conclude that under conservative assumptions 
only for values R < 0.8 a constraint on 7 can be derived 

Fleischer has argued that one can improve upon the above analysis by extracting 
some of the unknown hadronic parameters qc, s a , loq and rj from measurements 
of other decay processes [9]. The idea is to combine information on the ratio R 
with measurements of the direct CP asymmetries in the decays B° — > tt t K ± and 
B^ — > 7T ± K°, as well as of the ratio Rk defined in (2.29). One can then derive 
a bound on R that depends, besides the electroweak penguin parameters qc and 
u>c, only on a combination w = w(e a ,r)), which can be determined up to a two-fold 
ambiguity assuming SU(3) flavour symmetry. Besides the fact that this approach 
relies on SU(3) symmetry and involves significantly more experimental input than the 
original Fleischer-Mannel analysis, it does not allow one to eliminate the theoretical 
uncertainty related to the presence of electroweak penguin contributions. 



4. Hadronic uncertainties in the bound 

As a second application, we investigate the implications of recattering effects on 
the bound on COS7 derived from a measurement of the ratio R* defined in (1.2). 
In this case, the theoretical analysis is cleaner because there is model-independent 
information on the values of the hadronic parameters £3/2, q and u entering the 
parametrization of the isospin amplitude A 3 / 2 i n (2-5). The important point noted 
in [10] is that the decay amplitudes for B ± — » ir^K and 5 ± — > n^K^ differ only in 
this single isospin amplitude. Since the overall strength of A 3 / 2 is governed by the pa- 
rameter £3/2 and thus can be determined from experiment without much uncertainty, 
we have suggested to derive a bound on cos 7 without eliminating this parameter. In 
this respect, our strategy is different from the Fleischer-Mannel analysis. 

The exact theoretical expression for the inverse of the ratio R* is given by 

j COS0((5 E W - COS7) +£aCos(0-^)(l-<5 E wcos7) 

-ft* = 1 + ^£3/2 1 = 

yl - 2e a cost? cos 7 + e\ 

+ £3/2(1 - 25 EW cos7 + 5 EW ) , (4.1) 

where £3/2 has been defined in (2.6). Relevant for the bound on cos 7 is the maximal 
value R~ x can take for fixed 7. In [10], we have worked to linear order in the 
parameters £j, so that terms proportional to e a could be neglected. Here, we shall 
generalize the discussion and keep all terms exactly. Varying the strong-interaction 
phases and rj independently, we find that the maximum value of R" 1 is given by 

D-l ^ 1 , o- I^ew - COS7 ±£ a (l - 5 E wCOS7)| 2 n ox , X2 \ f An \ 

R* < l + 2£ 3 / 2 1 — + £ 3/2 (1-2(5ewcos7 + (5 ew ) , (4.2) 

yl =F2£ a cos7 + £2 
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where the upper (lower) signs apply if cos 7 < c (cos 7 > c ) with 

(1 + el) 5 E w x ( s 

c o = T—72T2 — - 5 ew • (4.3) 

' a °EW 

Keeping all terms in e 3 / 2 exactly, but working to linear order in e a , we find the 
simpler result 

K 1 < (l + £3/2 |5ew - cos 7 |) 2 + £3/2(^3/2 + 2e a ) sin 2 7 + 0(e 3/2 e 2 a ) . (4.4) 

The higher-order terms omitted here are of order 1% and thus negligible. The an- 
nihilation contribution e a enters this result in a very transparent way: increasing e a 
increases the maximal value of R^ 1 and therefore weakens the bound on cos 7. 

In [10], we have introduced the quantity A* by writing i?* = (1 — A*) 2 , so that 
A* = 1 — v^R* obeys the bound shown in (1.3). Note that to first order in e 3 / 2 the 
rescattering contributions proportional to e a do not enter. 6 Armed with the result 
(4.2), we can now derive the exact expression for the maximal value of the quantity 
A*, corresponding to the minimal value of R*. It is of advantage to consider the ratio 
A*/£ 3 /2, the bound for which is to first order independent of the parameter £ 3 / 2 . We 
recall that this ratio can be determined experimentally up to nonfactorizable SU(3)- 
breaking corrections. Its current value is A*/ 63/2 = 1.33 ± 0.78. 

In the left-hand plot in Figure 4, we show the maximal value for the ratio A*/£ 3 / 2 
for different values of the parameters e 3 / 2 and e a . The upper (red) and lower (blue) 
pairs of curves correspond to e 3 / 2 = 0.18 and 0.30, respectively, and span the allowed 
range of values for this parameter. For each pair, the dashed and solid lines corre- 
spond to e a = and 0.1, respectively. To saturate the bound (4.2) requires to have 
rj — <fi = 0° or 180°, in which case e a = 0.1 is a conservative upper limit (see Figure 1). 
The dotted curve shows for comparison the linearized result obtained by neglecting 
the higher-order terms in (1.3). The parameter <5ew = 0-64 is kept fixed in this plot. 
As expected, the bound on the ratio A*/£ 3 / 2 is only weakly dependent on the values 
of e 3 / 2 and e a . In particular, not much is lost by using the conservative value e a = 0.1. 
Note that for values A*/£ 3 / 2 > 0.8 the linear bound (1.3) is conservative, i.e., weaker 
than the exact bound, and even for smaller values of A*/e 3 / 2 the violations of this 
bound are rather small. Expanding the exact bound to next-to-leading order in f 3 / 2 , 
we obtain 

l +0(4 /2 ), (4.5) 



£3/2 



< \S E w - cos7| - e 3 / 2 



^3/2, 




showing that A^/e 3 / 2 > (1/2 + ^a/^3/2) 1 ^ 2 is a criterion for the validity of the lin- 
earized bound. This generalizes a condition derived, for the special case e a <^ e 3 / 2 , 
in [10]. 

6 Contrary to what has been claimed in [11], this does not mean that we were ignoring rescattering 
effects altogether. At linear order, these effects enter only through the strong-interaction phase 
difference which we kept arbitrary in deriving the bound on cos 7. 
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Figure 4: Left: Theoretical upper bound on the ratio versus I7I for various choices 

of £3/2 and e a , as explained in the text. The dotted line shows the linear bound derived 
in [10]. Right: Results for this ratio obtained for different choices of hadronic parameters. 
The curves show the theoretical bound for e a = 0.1 (solid) and e a = (dashed). The band 
shows the current experimental value with its la variation. 



To obtain a reliable bound on the weak phase 7, we must account for the the- 
oretical uncertainty in the value of the electroweak penguin parameter <5ew in the 
Standard Model, which is however straightforward to do by lowering (increasing) 
the value of this parameter used in calculating the right (left) branch of the curves 
defining the bound. The solid line in the right-hand plot in Figure 4 shows the most 
conservative bound obtained by using e a — 0.1 and varying the other two parameters 
in the ranges 0.18 < e 3 / 2 < 0.30 and 0.49 < <5ew < 0.79. The scatter plot shows the 
distribution of values of A*/ £3/2 obtained by scanning the strong-interaction param- 
eters over the same ranges as we did for the Fleischer-Mannel case in the previous 
section. The horizontal band shows the current central experimental value with its 
la variation. Unlike the Fleischer-Mannel bound, there is no violation of the bound 
(by construction), since all parameters are varried over conservative ranges. Indeed, 
for the points close to the right branch of the bound rj — = 0°, so that according to 
Figure 1 almost all of these points have e a < 0.03, which is smaller than the value we 
used to obtain the theoretical curve. The dashed curve shows the bound for e a = 0, 
which is seen not to be violated by any point. This shows that the rescattering effects 
parametrized by the quantity e a play a very minor role in the bound derived from the 
ratio R*. We conclude that, if the current experimental value is confirmed to within 
one standard deviation, i.e., if future measurements find that A*/£ 3 / 2 > 0.55, this 
would imply the bound I7I > 75°, which is very close to the value of 77° obtained in 
[10]. 
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25. 50. 75. 100. 125. 150. 175. 25. 50. 75. 100. 125. 150. 175. 

gamma gamma 

Figure 5: Left: Results for the ratio R* obtained for different choices of hadronic param- 
eters. The curves show the theoretical bound for e a = 0.1 (solid) and e a = (dashed). 
The band shows the current experimental value of R* with its la variation. Right: Results 
for the ratio R* for different values of the strong- interaction phase \4>\, as indicated by the 
numbers. 



Given that the experimental determination of the parameter £3/2 is limited by 
unknown nonfactorizable SU(3)-breaking corrections, one may want to be more con- 
servative and derive a bound directly from the measured ratio R* rather than the 
ratio In the left-hand plot in Figure 5, we show the same distribution as 

in the right-hand plot in Figure 4, but now for the ratio R*. The resulting bound 
on 7 is slightly weaker, because now there is a stronger dependence on the value of 
£3/2, which we vary as previously between 0.18 and 0.30. If the current value of R* 
is confimed to within one standard deviation, i.e., if future measurements find that 
R* < 0.71, this would imply the bound I7I > 72°. 

Besides providing interesting information on 7, a measurement of R* or A*/e* 
can yield information about the strong-interaction phase <fi. In the right plot in 
Figure 5, we show the distribution of points obtained for fixed values of the strong- 
interaction phase |0| between 0° and 180° in steps of 30°. For simplicity, the param- 
eters £3/2 = 0.24 and <5ew = 0.64 are kept fixed in this plot, while all other hadronic 
parameters are scanned over the realistic parameter set. We observe that, indepen- 
dently of 7, a value R* < 0.8 requires that \<p\ < 90°. This conclusion remains true 
if the parameters £3/2 and <5ew are varied over their allowed ranges. We shall study 
the correlation between the weak phase 7 and the strong phase in more detail in 
Section 6. 

Finally, we emphasize that a future, precise measurement of the ratio R* may 
also yield a surprise and indicate physics beyond the Standard Model. The global 
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analysis of the unitarity triangle requires that |7| < 105° [32], for which the lowest 
possible value of R* in the Standard Model is about 0.55. If the experimental value 
would turn out to be less than that, this would be strong evidence for New Physics. 
In particular, in many extensions of the Standard Model there would be additional 
contributions to the electroweak penguin parameter <5ew arising, e.g., from penguin 
and box diagrams containing new charged Higgs bosons. This could explain a larger 
value of R*. Indeed, from (4.4) we can derive the bound 

^R- 1 - £3/2(^3/2 + 2£ a ) sin 2 7 max - 1 
Oew > z 

£3/2 

+ cos7 max , (4.6) 

where 7 max is the maximal value allowed by the global analysis (assuming that 7 max > 
arccos(co) ~ 50°). In Figure 6, we show this bound for the current value 7 max = 105° 
and three different values of £3/2 as well as two different values of e a . The gray band 
shows the allowed range for <5ew hi the Standard Model. In the hypothetical situation 
where the current central values i?* = 0.47 and £3/2 = 0.24 would be confirmed by 
more precise measurements, we would conclude that the value of 5ew is at least twice 
as large as predicted by the Standard Model. 



5. Prospects for direct CP asymmetries and prediction for 
the 

B° — > 7T°K° branching ratio 

5.1 Decays of charged B mesons 

We will now analyse the potential of the various B — > irK decay modes for show- 
ing large direct CP violation, starting with the decays of charged B mesons. The 
smallness of the rescattering effects parametrized by e a (see Figure 1) combined with 
the simplicity of the isospin amplitude A3/2 (see Section 2.3) make these processes 
particularly clean from a theoretical point of view. 

Explicit expressions for the CP asymmetries in the various decays can be de- 
rived in a straightforward way starting from the isospin decomposition in (2.2) and 
inserting the parametrizations for the isospin amplitudes derived in Section 2. The 
result for the CP asymmetry in the decays B ± — > ^K has already been presented 
in (2.27). The corresponding expression for the decays B ± — > ir ^ reads 

a t o • d £3/2sin0 + £ a sin77-5 3/ 2£a<5Ewsin(0-r7) 

A CP 7T K + ) = 2sm7i?, —!■ '- — , 5.1 

1 — 2e a COS T) cos 7 + £q 

where the theoretical expression for R* is given in (4.1), and we have not replaced 
£3/2 in terms of e 3 / 2 . Neglecting terms of order e a and working to first order in 
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£3/2, we find the estimate A C p(n K + ) ~ 2e 3 / 2 sin 7 sin <fi ps 0.5 sin 7 sin 0, indicating 
that potentially there could be a very large CP asymmetry in this decay (note that 
sin 7 > 0.73 is required by the global analysis of the unitarity triangle). 
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Figure 7: Results for the direct CP asymmetries Acp(tt + K°) (left) and Acp(tt°K + ) 
(right) versus |7|. 

In Figure 7, we show the results for the two direct CP asymmetries in (2.27) and 
(5.1), both for the realistic and for the conservative parameter sets. These results 
confirm the general observations made above. For the realistic parameter set, and 
with 7 between 47° and 105° as indicated by the global analysis of the unitarity 
triangle [32], we find CP asymmetries of up to 15% in 5 ± — > ^K decays, and of 
up to 50% in B^ — > tt°K ± decays. Of course, to have large asymmetries requires 
that the sines of the strong-interaction phases 77 and <p are not small. However, this is 
not unlikely to happen. According to the left-hand plot in Figure 1, the phase r) can 
take any value, and the phase could quite conceivably be large due to the different 
decay mechanisms of tree- and penguin-initiated processes. We stress that there is no 
strong correlation between the CP asymmetries in the two decay processes, because 
as shown in Figure 1 there is no such correlation between the strong-interaction 
phases r\ and 0. 

5.2 Decays of neutral B mesons 

Because of their dependence on the hadronic parameters £t, 1c and ojc entering 
through the sum A 1 / 2 + A 3 / 2 of isospin amplitudes, the theoretical analysis of neutral 
B — > txK decays is affected by larger hadronic uncertainties than that of the decays 
of charged B mesons. Nevertheless, some interesting predictions regarding neutral 
B decays can be made and tested experimentally. 
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The expression for the direct CP asym- 
metry in the decays B° — > 7r Zf K ± 



is 



2 sin 7 St (sin — e T q c sin u c ) \ § 




i2 1 
(5.2) 

where = 0t — 0p- This result reduces to 
(5.1) under the replacements qc — > <5ew, 

— > 0, i? — > R" 1 , and — > £ 3 / 2 . The 
corresponding expression for the direct 
CP asymmetry in the decays B° — > tc°K° 
and _B° — > 7T°K° is more complicated and 
will not be presented here. Below, we 
shall derive an exact relation between the 
various asymmetries, which can be used 
to compute Acp(n°K°). 

Gronau and Rosner have emphasized 
that one expects Acp(tt^K + 

and that one could thus combine the data spectively. The solid and dashed lines refer 
samples for these decays to enhance the to e a = 0.1 and 0. The band shows the pre- 
statistical significance of an early signal diction for Sew in the Standard Model, 
of direct CP violation [36] . We can easily 

understand the argument behind this observation using our results. Neglecting the 
small rescattering contributions proportional to e a for simplicity, we find 



Figure 6: Lower bound on the parameter 
Sew as a function of R*. The upper (red), 
middle (black) and lower (blue) curves cor- 
A CP (tt K + Respond to e 3/2 = 0.18, 0.24 and 0.30, re- 



A C p(ti K + ) 1 e r (sin0 — e r gc sin c^c) 1 Et 
A CP (ir K+) ~ R*R £ 3/ / 2 sin ~~ R*R e 3 / 2 



(5.3) 



In the last step, we have used that the electroweak penguin contribution is very 
small because it is suppressed by an additional factor of et, and that the strong- 
interaction phases <ft and <fi are strongly correlated, as follows from the right-hand 
plot in Figure 1. Numerically, the right-hand side turns out to be close to 1 for 
most of parameter space. This is evident from the left-hand plot in Figure 8, which 
confirms that there is indeed a very strong correlation between the CP asymmetries 



in the decays B c 



tt^K* and B ± 



7T 



in agreement with the argument 



given in [36]. Combining the data samples for these decays collected by the CLEO 
experiment, one may have a chance for observing a statistically significant signal for 
the first direct CP asymmetry in B decays before the operation of the asymmetric 
B factories. 



The decays B c 



7c°K° and B° 



7T 



'K have not yet been observed experi- 



mentally, but the CLEO Collaboration has presented an upper bound on their CP- 
averaged branching ratio of 4.1 x 10~ 5 [1]. In analogy with (1.2), we define the 
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Figure 8: Correlation between the direct CP asymmetries ACPI = Acp(tt°K + ), ACP2 = 
A CP (tt-K+), and ACP3 = A C p(vr ^ ). 

ratios 

r(B+) 2[Br(ff° -> n°K°) + Bi(B° -> 7r°i?°)] 
~ t(B°) Br(B+ TT+iT ) + Br(£- -> 7^7?°) ' 

2[Br(£?° -> 7r°K°) + Br(E° -> 7r°i?°)] _ # 
°* ~ Br(5° -> tt-^+) + Br(5° -> tt+^-) ~ i? ' ( ' 

Using our parametrizations for the different isospin amplitudes, we find that the 
ratios R, i?* and Ro obey the relations 

Ro-R + R; 1 -1 = A U R -RR* = A 2 + O(e 3 i ), (5.5) 

where 

Ai = 2e\ j2 (1 - 2fe w cos 7 + <J| W ) - 2e 3/2 e T (f - fe w cos 7) cos(0 T - 3/2 ) 

- 2e 3/2 e T qc (few - cos 7) cos(0 T - 03/2 + v c ) , 
A 2 = A x — 4^3/2 (fe w — COS7) 2 cos 2 

+ 4e 3 / 2 £t (few _ cos 7) cos0 qc cos(0 + uc) — cos 7 cos0 , (5.6) 

and £t is defined in analogy with e 3 / 2 in (2.6), so that Et/st = ^3/2/^3/2- The first 
relation in (5.5) generalizes a sum rule derived by Lipkin, who neglected the terms 
of 0{ef) on the right-hand side as well as electroweak penguin contributions [37]. 
The second relation is new. It follows from the fact that i?o* = R* + 0(e 2 ), which is 
evident since the pairs of decay amplitudes entering the definition of the two ratios 
differ only in the isospin amplitude A 3 / 2 . 



25 






25. 50. 75. 100. 125. 150. 175. 0.4 0.6 0.8 1. 1.2 1.4 1.6 

gamma Rstar 

Figure 9: Ratio Rq (left), and correlation between the ratios R* and i?o* (right). 



The left-hand plot in Figure 9 shows the results for the ratio Rq versus |7|. 
The dependence of this ratio on the weak phase turns out to be much weaker than 
in the case of the ratios R and R*. For the realistic parameter set we find that 
0.7 < Rq < 1.0 for most choices of strong- interaction parameters. Combining this 
with the current value of the — > n^K branching ratio, we obtain values between 
(0.47±0. 18) x 10~ 5 and (0.67±0.26) x 10~ 5 for the CP-averaged B° -> ti K branching 
ratio. The right-hand plot in Figure 9 shows the strong correlation between the ratios 
R* and i? * — Ro/R, which holds with a remarkable accuracy over all of parameter 
space. 

In Figure 10, we show the estimates of Rq obtained by neglecting the terms of 
0{ef) and higher in the two sum rules in (5.5). Using the present data for the various 
branching ratios yields to the estimates Rq = (—0.1 ± 0.9) from the first and Rq = 
(0.5 ±0.2) from the second sum rule. Both results are consistent with the theoretical 
expectations for Rq exhibited in the left-hand plot in Figure 9; however, the second 
estimate has a much smaller experimental error and, according to Figure 10, it is 
likely to have a higher theoretical accuracy. We can rewrite this estimate as 

i[B,(B° - WBr( B° - ■»*«)] * ^ ■ M 

where the branching ratios on the right-hand side are averaged over CP-conjugate 
modes. With current data, this relation yields the value (0.33 ± 0.18) x 10" 5 . Com- 
bining the three estimates for the CP-averaged B° — > n K° branching ratio presented 
above we arrive at the value (0.5 ± 0.2) x 10~ 5 , which is about a factor of 3 smaller 
than the other three B — > nK branching ratios quoted in (1.1). 

We now turn to the study of the direct CP asymmetry in the decays B° — > 7i°K° 
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Figure 10: Esimates of Ro obtained using the first (left) and second (right) sum rule in 
(5.5). 

and B° — > n°K°. Using our general parametrizations, we find the sum rule 

A CP (ir + K°) - R; 1 A CP (n°K + ) + R A^- K+) - R A CP (ir°K°) 

= 2 sin 7 e 3/2 e T [# E w sin(0 T - 3/2 ) - q c sin(0 T - 3/2 + uj c )\ • (5.8) 

By scanning all strong-interaction parameters, we find that for the realistic (conser- 
vative) parameter set the right-hand side takes values of less that 4% (7%) times sin 7 
in magnitude. Neglecting these small terms, and using the approximate equality of 
the CP asymmetries in B ± — > tt°K ± and B° — > n T K ± decays as well as the second 
relation in (5.5), we obtain 

A 1 - RR * A fT Q^ 1 ^Cp(vt + K°) 

A C p(7t A ) ~ A cp (ti K ) + — . (5.9) 

it JX^ it rt* 

The first term is negative for most choices of parameters and would dominate if the 
CP aymmetry in B ± — > 7r°i^ ± decays would turn out to be large. We therefore 
expect a weak anticorrelation between Acp(ir°K ) and A C p(7r + i^ ), which is indeed 
exhibited in the right-hand plot in Figure 9. 

For completeness, we note that in the decays B°, B° — > n°Ks one can also study 
mixing-induced CP violation, as has been emphasized recently in [11]. Because of 
the large hadronic uncertainties inherent in the calculation of this effect, we do not 
study this possibility further. 

6. Determination of 7 from B ± — > nK, 7T7T decays 

Ultimately, one would like not only to derive bounds on the weak phase 7, but to 
measure this parameter from a study of CP violation in B — > tcK decays. However, as 
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we have pointed out in Section 2, this is not a trivial undertaking because even perfect 
measurements of all eight B — > ixK branching ratios would not suffice to eliminate 
all hadronic parameters entering the parametrization of the decay amplitudes. 

Because of their theoretical cleanness, the decays of charged B mesons are best 
suited for a measurement of 7. In [14], we have described a strategy for achieving this 
goal, which relies on the measurements of the CP-averaged branching ratios for the 
decays B ± — > tt^K and B^ — > tt^tt , as well as of the individual branching ratios 
for the decays B + — > tt°K + and B~ — > tt°K~, i.e., the direct CP asymmetry in 
this channel. This method is a generalization of the Gronau-Rosner-London (GRL) 
approach for extracting 7 [38]. It includes the contributions of electroweak penguin 
operators, which had previously been argued to spoil the GRL method [29, 39]. 

The strategy proposed in [14] relies on the dynamical assumption that there 
is no CP-violating contribution to the B ± — > n^K decay amplitudes, which is 
equivalent to saying that the rescattering effects parametrized by the quantity e a in 
(2.4) are negligibly small. It is evident from the left-hand plot in Figure 1 that this 
assumption is indeed justified in a large region of parameter space. Here, we will 
refine the approach and investigate the theoretical uncertainty resulting from e a 7^ 0. 
As a side product, we will show how nontrivial information on the strong-interaction 
phase difference = 03/2 — <pp can be obtained along with information on 7. 

To this end, we consider in addition to the ratio the CP-violating observable 

t A CP (ir°K + ) + _ sin - e a 5 EW sin(0 - rj) 

A= Acp{n + K v ) = 2 sin 7e 3/2 , =■ (6.1) 

-"-* J 1 — 2e a cos i] cos 7 + el 

The purpose of subtracting the CP asymmetry in the decays B ± -> ^K is to elim- 
inate the contribution of 0(e a ) in the expression for Acp(tt°K + ) given in (5.1). A 
measurement of this asymmetry is the new ingredient in our approach with respect 
to that in [14]. With the definition of A as given above, the rescattering effects 
parametrized by e a are suppressed by an additional factor of e 3 / 2 and are thus ex- 
pected to be very small. As shown in Section 4, the same is true for the ratio R*. 
Explicitly, we have 

R* 1 = 1 + 2e 3/2 cos0 (5 E w - C0S7) +£3/2 (1 - 2<5 E wcos7 + <5| w ) + 0(^3/2^) , 

A = 2 sin 7^3/2 sin + 0(e 3/2 e a ) ■ (6.2) 

These equations define contours in the (7, 0) plane. When higher-order terms 
are kept, these contours become narrow bands, the precise shape of which depends 
on the values of the parameters e 3 / 2 and 5ew- 111 the limit e a = the procedure 
described here is mathematically equivalent to the construction proposed in [14]. 
There, the errors on cos 7 resulting from the variation of the input parameters have 
been discussed in detail. For a typical example, where 7 = 76° and = 20°, we 
found that the uncertainties resulting from a 15% variation of e 3 / 2 and <5ew are 
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cos 7 = 0.24 ± 0.09 ± 0.09, correspondig to errors of ±5° each on the extracted value 
of 7. 

Our focus here is to evaluate the additional uncertainty resulting from the rescat- 
tering effects parametrized by e a and 77. For given values of e 3 / 2 , few, £a, ?7, an d 7, 
the exact results for R* in (4.1) and A in (6.1) can be brought into the generic form 
A cos + B sin = C, where in the case of 

A few - COS7 + £ a cosr? (1 -few cos 7) 

A = Z€ 3 / 2 



B = 2e a 



yl - 2e a cos 77 cos 7 + el 

e a smr] (1 - fe w cos 7) 
3 / 2 n ~ ? ' 



whereas for A 



1 1 — 2e a cos 77 cos 7 + e\ 
C = R; 1 -1- e\ j2 (1 - 2fe w cos 7 + <S| W ) , (6.3) 



n _ e a <5 EW sin 7/ 



-B = 2£ 3 /2 



\Jl — 2e a cos 77 cos 7 + 

1 - e a few cost; 
^1 - 2e a cos 77 cos 7 + 



C = — . (6.4) 
sin 7 

The two solutions for cos0 are given by 

, AC±BVA2 + B*-C2 
cos0= A2 + B2 . (6.5) 

The physical solutions must be such that cos0 is real and its magnitude less than 1. 

In Figure 12, we show the resulting contour bands obtained by keeping e 3 / 2 = 
0.24 and few = 0.64 fixed to their central values, while the rescattering parameters 
are scanned over the ranges < e a < 0.08 and —180° < 77 < 180°. Assuming that 
sin 7 > as suggested by the global analysis of the unitarity triangle, the sign of A 
determines the sign of sin0. In the plot, we assume without loss of generality that 
0° < <p < 180°. For instance, if i?* = 0.7 and A = 0.2, then the two solutions are 
(7, 0) w (98°, 25°) and (7, 0) w (153°, 67°), only the first of which is allowed by the 
upper bound 7 < 105° following from the global analysis of the unitarity triangle 
[32]. It is evident that the contours are rather insensitive to the rescattering effects 
parametrized by e a and 77. The error on 7 due to these effects is about ±5°, which 
is similar to the errors resulting from the theoretical uncertainties in the parameters 
e 3 /2 and few- The combined theoretical uncertainty is of order ±10° on the extracted 
value of 7. 
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To summarize, the strategy for deter- 
mining 7 would be as follows: From mea- 
surements of the CP-averaged branching 
ratio for the decays B — > n^n , — > 
vr ± K° and B ± -> vr ^, the ratio and 
the parameter e 3 / 2 are determined using 
(1.2) and (1.4), respectively. Next, from 
measurements of the rate asymmetries in 
the decays B ± n ± K° and B ± n ^ 
the quantity A is determined. From the 
contour plots for the quantities i?* and 
A the phases 7 and <fi can then be ex- 
tracted up to discrete ambiguities. In 
this determination one must account for 
theoretical uncertainties in the values of 
the parameters e 3 / 2 and 5ew ; as well as 
for rescattering effects parametrized by e a 
and rj. Quantitative estimates for these 
uncertainties have been given above. 

7. Conclusions 
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Figure 11: Contour plots for the quanti- 
ties i?* (red "hyperbolas") and A (blue "cir- 
cles"). The scatter plots show the results 
including rescattering effects, while the lines 
refer to e a = 0. The solid curves correspond 
to the contours for = 0.7 and A = 0.2, 
the dashed ones to R* = 0.9 and A = 0.4. 



We have presented a model-independent, global analysis of the rates and direct CP 
asymmetries for the rare two-body decays B — > nK. The theoretical description 
exploits the flavour symmetries of the strong interactions and the structure of the 
low-energy effective weak Hamiltonian. Isospin symmetry is used to introduce a 
minimal set of three isospin amplitudes. The explicit form of the effective weak 
Hamiltonian in the Standard Model is used to simplify the isovector part of the 
interaction. Both the numerical smallness of certain Wilson coefficient functions 
and the Dirac and colour structure of the local operators are relevant in this con- 
text. Finally, the £7-spin subgroup of flavour SU(3) symmetry is used to simplify the 
structure of the isospin amplitude A 3 / 2 referring to the decay B — > {tiK) i=3 / 2 . In 
the limit of exact [/-spin symmetry, two of the four parameters describing this am- 
plitude (the relative magnitude and strong-interaction phase of electroweak penguin 
and tree contributions) can be calculated theoretically, and one additional parame- 
ter (the overall strength of the amplitude) can be determined experimentally from a 
measurement of the CP-averaged branching ratio for B ± — > ty^tt decays. What re- 
mains is a single unknown strong-interaction phase. The SU(3)-breaking corrections 
to these results can be calculated in the generalized factorization approximation, so 
that theoretical limitations enter only at the level of nonfactorizable SU(3)-breaking 
effects. However, since we make use of SU(3) symmetry only to derive relations for 
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amplitudes referring to isospin eigenstates, we do not expect gross failures of the 
generalized factorization hypothesis. We stress that the theoretical simplifications 
used in our analysis are the only ones rooted on first principles of QCD. Any further 
simplification would have to rest on model-dependent dynamical assumptions, such 
as the smallness of certain flavour topologies with respect to others. 

We have introduced a general parametrization of the decay amplitudes, which 
makes maximal use of these theoretical constraints but is otherwise completely gen- 
eral. In particular, no assumption is made about strong-interaction phases. With the 
help of this parametrization, we have performed a global analysis of the branching 
ratios and direct CP asymmetries in the various B — > irK decay modes, with par- 
ticular emphasis on the impact of hadronic uncertainties on methods to learn about 
the weak phase 7 = &rg(V* b ) of the unitarity triangle. The main phenomenological 
implications of our results can be summarized as follows: 

• There can be substantial corrections to the Fleischer-Mannel bound on 7 from 
enhanced electroweak penguin contributions, which can arise in the case of 
a large strong-interaction phase difference between I — \ and I — \ isospin 
amplitudes. Whereas these corrections stay small (but not negligible) if one 
restricts this phase difference to be less than 45°, there can be large violations 
of the bound if the phase difference is allowed to be as large as 90°. 

• On the contrary, rescattering effects play a very minor role in the bound on 7 
derived from a measurement of the ratio of CP-averaged B^ 1 — > 7iK branch- 
ing ratios. They can be included exactly in the bound and enter through a 
parameter e a , whose value is less than 0.1 even under very conservative con- 
ditions. Including these effects weakens the bounds on 7 by less than 5°. We 
have generalized the result of our previous work [10], where we derived a bound 
on cos 7 to linear order in an expansion in the small quantity £ 3 / 2 . Here we 
refrain from making such an approximation; however, we confirm our previous 
claim that to make such an expansion is justified (i.e., it yields a conservative 
bound) provided that the current experimental value of R* does not change 
by more than one standard deviation. The main result of our analysis is given 
in (4.2), which shows the exact result for the maximum value of the ratio R* 
as a function of the parameters <5ew> £3/2, and e a . The first parameter de- 
scribes electroweak penguin contributions and can be calculated theoretically. 
The second parameter can be determined experimentally from the CP-averaged 
branching ratios for the decays B ± — > 7r ± 7r° and B ± — > ^K . We stress that 
the definition of e 3 / 2 is such that it includes exactly possible rescattering contri- 
butions to the B ± — > 7T^~K® decay amplitudes. The third parameter describes a 
certain class of rescattering effects and can be constrained experimentally once 
the CP-averaged B^ — > K ± K° branching ratio has been measured. However, 
we have shown that under rather conservative assumptions e a < 0.1. 
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• The calculable dependence of the B ± — > 7T.K" decay amplitudes on the elec- 
troweak penguin contribution 5 EW offers a window to New Physics. In many 
generic extensions of the Standard Model such as multi-Higgs models, we ex- 
pect deviations from the value <5ew — 0-64 ± 0.15 predicted by the Standard 
Model. We have derived a lower bound on <5ew as a function of the value of 
the ratio and the maximum value for 7 allowed by the global analysis of the 
unitarity triangle. If it would turn out that this value exceeds the Standard 
Model prediction by a significant amount, this would be strong evidence for 
New Physics. In particular, we note that if the current central value i?* = 0.47 
would be confirmed, the value of 5ew would have to be at least twice its stan- 
dard value. 

• We have studied in detail the potential of the various B — > nK decay modes 
for showing large direct CP violation and investigated the correlations between 
the various asymmetries. Although in general the theoretical predictions suffer 
from the fact that an overall strong-interaction phase difference is unknown, we 
conclude that there is a fair chance for observing large direct CP asymmetries 
in at least some of the decay channels. More specifically, we find that the direct 
CP asymmetries in the decays —> and B° — > tt^ 1 K ± are almost fully 
correlated and can be up to 50% in magnitude for realistic parameter choices. 
The direct CP asymmetry in the decays B° — > 7r°K° and B° — > ir°K° tends 
to be smaller by about a factor of 2 and anticorrelated in sign. Finally, the 
asymmetry in the decays B^ — > n^K is smaller and uncorrelated with the 
other asymmetries. For realistic parameter choices, we expect values of up to 
15% for this asymmetry. 

• We have derived sum rules for the branching ratio and direct CP asymmetry 
in the decays B° — > 7i°K° and B° — > tt°K . A rather clean prediction for the 
CP-averaged branching ratio for these decays in given in (5.7). We expect a 
value of (0.5 ± 0.2) x 10~ 5 for this branching ratio, which is about a factor of 
3 less than the other B — > nK branching ratios. 

• Finally, we have presented a method for determining the weak phase 7 along 
with the strong-interaction phase difference from measurements of B ± — > %K, 
-kit branching ratios, all of which are of order 10~ 5 . This method generalizes 
an approach proposed in [14] to include rescattering corrections to the B^ — > 
tt^K decay amplitudes. We find that the uncertainty due to rescattering 
effects is about ±5° on the extracted value of 7, which is similar to the errors 
resulting from the theoretical uncertainties in the parameters £3/2 and Sew- 
The combined theoretical uncertainty in our method is of order ±10°. 

A global analysis of branching ratios and direct CP asymmetries in rare two-body 
decays of B mesons can yield interesting information about fundamental parameters 
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of the flavour sector of the Standard Model, and at the same time provides a window 
to New Physics. Such an analysis should therefore be a central focus of the physics 
program of the B factories, which in many respects is complementary to the time- 
dependent studies of CP violation in neutral B decays into CP eigenstates. 
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